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at rest coalesce to form a stable bound soliton molecule due to dipolar attraction. 
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A bright soliton is a self-reinforcing solitary wave that maintains its shape, while 
traveling at constant speed, due to a cancellation of nonlinear attraction and dispersive 
effects. Experimentally, bright matter- wave solitons and soliton trains were created in a 
Bose-Einstein condensate (BEC) of 7 Li [TJ [2] and 85 Rb atoms [3] by turning the atomic 
interaction attractive from repulsive using a Feshbach resonance (FR) [4J and releasing 
the BEC in an axially free or an expulsive trap. 

Lately, a BEC of 52 Cr atoms with a large long-range dipolar interaction has been 
made (5J [6]. This allows to study the dipolar BEC (DBEC) with variable short-range 
interaction j5j [6] using a FR [4] . The DBEC has many distinct features [5J U\ El E] . 
The stability of a DBEC depends not only on the scattering length, but also on 
the trap geometry [9j. A disk-shaped trap gives a repulsive dipolar interaction 
and the DBEC is more stable. In contrast, a cigar-shaped trap yields an attractive 
dipolar interaction and hence may favor a collapse [5j [9j [TO]. Collapse instability 
of solitons in the nonpolynomial Schrodinger equation with dipole-dipole interactions 
was studied by Gligoric et al [H]. The controllable short-range interaction with the 
exotic dipolar interaction makes the DBEC an interesting system for soliton generation 
in the cigar configuration due to added dipolar attraction and a challenging system for 
theoretical investigation [12] to study the interplay between the dipolar and short-range 
interactions. 

We suggest cigar-shaped bright and vortex solitons [13j for repulsive short-range 
interaction in a DBEC free to move along the axial direction and trapped in the radial 
direction, using the three-dimensional (3D) Gross-Pitaevskii (GP) equation, and study 
their statics and dynamics. Vortex-bright solitons (called vortex solitons hereafter) 
are solitons of (quantized) unit angular momentum. In the usual (non-dipolar) BEC, 
solitons with intrinsic vorticity trapped in the cigar-shaped potential, were considered 
by Salasnich et al [14J. Vortex states of higher angular momentum (Z > 1) are unstable 
and decay into multiple vortices of angular momentum ft [15j. For a conventional 
BEC without dipole moment, solitons appear only for attractive atomic interaction. 
Phase plots showing the stability of the solitons for different dipolar and short-range 
interactions are obtained. We study breathing oscillation and collision of two bright 
and two vortex solitons as well as the statics and dynamics of the bright solitons using 
a Gaussian variational approach. Collision of two bright and vortex solitons in non- 
dipolar BEC was studied in [16J. Where applicable, the variational results are found 
to be in good agreement with numerical results. By reversing the sign of the dipolar 
interaction, Pedri and Santos [FT] studied stable two-dimensional (2D) bright solitons in 
a disk-shaped DBEC free to move in the radial plane and trapped in the axial direction. 
These solitons were later generalized to 2D vortex solitons [T8] . There has also been 
study of anisotropic 2D solitons not requiring the sign-reversal of the dipolar interaction 
[T9] . However, special conditions are needed to prepare these 2D solitons, whereas the 
present one-dimensional (ID) solitons can be formed under very general conditions and 
hence, are of greater experimental interest. Solitons in the Tonks-Girardeau gas with 
dipolar interaction was studied by Baizakov et al [20j in a ID model. 
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Cuevas et al [21] studied solitons in cigar-shaped DBEC using a reduced quasi-lD 
toy model [22J. As there is no collapse in ID (models with cubic nonlinearity as in 
this case), such study predicts soliton for attractive interaction for an arbitrarily large 
number of atoms. But non-dipolar BEC bright solitons [23] and vortex solitons [21] in 
3D are stable upto a critical number of atoms, beyond which the system collapses. The 
same is found to be true for a DBEC; viz. the critical numbers reported in figure [j] (a) 
DBEC soliton. 

We study the bright solitons in a DBEC of TV atoms, each of mass m, using the 
GP equation: [5j 



.d<j>(r,t) 



dt 

+ 
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N J U dd (r-v')\<f>(r',t)\ 2 dT' 

with the dipolar interaction [/^(R) = 3a^(l — 3cos 2 #)/i? 3 , R = r — r' and 
normalization / (j)(r) 2 dr = 1. Here, a is the atomic scattering length, 6 the angle 
between R and the polarization direction z. The length add = Pop 2 m/ (12tt^ 2 ) measures 
the strength of dipolar interaction and its experimental value for 52 Cr is 15a [5], with 
a the Bohr radius, p the (magnetic) dipole moment of a single atom, and p the 
permeability of free space. The trap p 2 /2 acts in the radial p = (x, y) direction only. 
To obtain a quantized vortex of angular momentum Ih, I = 1, around axial z axis, one 
has to introduce a phase (equal to the azimuthal angle) in the wave function [25] . This 
procedure introduces a centrifugal term l 2 /2p 2 in the GP equation for a vortex and we 
adopt this method to study vortex solitons. In ([j]) we use oscillator unit of transverse 
trap: unit of time is inverse angular frequency uj± taken as 1 ms, unit of length is 



oscillator length a± = yh/muj^ taken as 1 pm. 

Lagrangian density of ([j]) for bright soliton is [13] 

' ' 1 * - + \\ V0| 2 + ^|0| 2 + 2iraN\<p\ 4 



2' r| 2 

+ yl0| 2 /Mr-r')|0(r')|W. (2) 

For a variational study we use the Gaussian ansatz [26j: 0(r, t) = exp(—p 2 /2w 2 — z 2 /2wl 
+iap 2 + if3z 2 )/(w p y^w z 7r 3 / 4 ) where w p and w z are time-dependent widths and a and /3 
are time-dependent phases. The effective Lagrangian L (per particle) is 



L = J Cdr = (w 2 p a + ^1/3 + 2w 2 a 2 + w 2 z f3 2 



with £ dip = N[a - a dd f(K)]/(V2^w 2 p w z )J(K) = [1 + 2k 2 - 3n 2 d( K )}/(l - n 2 ),d{n) = 
(atanhVl — ^ 2 )/Vl — k 2 , k = w p /w z . In the cigar shape, the dipolar interaction 
becomes attractive and contributes negatively in Sd\ p - The Euler- Lagrange equations 
for parameters w p ,w z ,a, f3 can be used to obtain the following equations for the widths 

1 1 N In 

w p + w p = — + —= — — [2a - a dd g(K)\ , (4) 
w 6 p V2irw*w z 




Figure 1. (a) iV cr it versus a plot for add — 15ao, lOOao- (b) rms sizes (p), (z) and 
chemical potential \i versus a for N = 1000 and add — lOOao- The two arrows denote 
the region of stability of bright solitons. (variational: v, numerical: n, bright: bri, 
vortex: vor.) 
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w z = — + —f=^-^ [a - a dd h(K)\ , (5) 
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with^(«) = [2-7K 2 -^ + 9^d{K)]/{l-K 2 ) 2 ,h{K) = [l + l()K 2 -2^-§K 2 d{K)]{l-K 2 ) 2 . 
Equations Q and ^ determine the dynamics of bright solitons. The stationary widths 
are obtained from these equations by setting w p = w z = 0. The stationary chemical 
potential (/x) is given by \i — l/2w 2 + l/4w^ + w 2 /2 + 2£ dip . 

We perform a 3D numerical simulation in Cartesian x,y, z variables employing 
imaginary- and real-time propagation with Crank-Nicolson method using Fortran 
programs provided in Ref. [27]. The dipolar interaction is evaluated by the usual 
fast Fourier transformation [7j. 

For bright solitons the system must be attractive. In the absence of dipole moment 
ifldd = 0), attraction corresponds to a < 0. For a d d > 0, the dipolar interaction in ([!]) 
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Figure 2. A sectional view of the 3D contour plot of density for a (a) bright 
(add = 15ao and a = 0.5 nm) and (b) vortex (add = lOOao, a = 4 nm, I = 1) 
solitons of N — 1000 atoms. The density on the contour is 0.001. 

contributes attractively in the cigar shape, which can compensate for the atomic short- 
range repulsion, thus forming bright solitons for a > 0. From ([5]), the onset of attraction 
for a bright soliton is at [a—dddhfo)] = while one must have TV oo. Then from Q, as 
TV — > oc, one must have [2a — a^g(ft)] = 0. These two conditions can be simultaneously 
satisfied as k for a weakly bound soliton of infinite extension along axial direction 
while h(n = 0) = g{n = 0)/2 = 1. Consequently, bright solitons are allowed for 
— oc < a < add and for TV below a critical number A^ crit , as ^ and ^ only permit 
solution for TV < A^ crit . The system collapses for TV > A^ crit . These results are illustrated 
in figure [l] (a), where we show A^ cri t versus a for add — 15a ( 52 Cr atom), and 100a , 
for bright and vortex solitons from numerical and variational calculation. The value 
add = 100a is considered to simulate a possible DBEC of dipolar atoms (or molecules) 
of larger dipole moment, such as 164 Dy 66 , or 166 Er 68 , which are being considered for 
BEC experiment [281 [29]. Preliminary study [29] indicates that the dysprosium atom 
has a dipole moment of about nine times larger than that of 52 Cr, which justifies the 
use of add = 100a . The attractive region of collapse (TV > N crit ,a < a^), stable 
soliton formation (N < A^ crit ,a < add) and the repulsive region of expansion (a > add) 
are clearly shown in figure [l] (a). Due to the centrifugal repulsion, the vortex soliton 
can accommodate a larger number of atoms [24]. In our numerical simulation we shall 
consider DBEC with add > a > — domain of soliton controlled solely by dipolar 
interaction and inaccessible for a normal BEC with add = 0. In figure [I] (b) we plot 
results for root-mean-square (rms) sizes (p), (z) and chemical potential \i of bright and 
vortex solitons. The variational results for bright solitons are in good agreement with 
the numerical results. From figures |l|(a) and (b) we find that the vortex soliton is stable 
in a wider domain of parameter space, due to the vortex core which keeps the atoms 
apart to avoid the collapse. Consequently, the vortex soliton is larger in size than the 
bright soliton with a larger chemical potential. 

To illustrate the shape of the quasi-lD soliton, in figure [2] (a) we show a section of 
the 3D contour plot of density of a bright soliton of 1000 atoms for a = 0.5 nm, and 
add — 15ao = 0.7938 nm. A similar plot for the vortex soliton for TV = 1000, a = 4 nm, 
and add = 100a = 5.2917 nm is shown in figure [2] (b) which clearly exhibits the vortex 
core along the z axis. 
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Next we study by numerical (num) and variational (var) approaches the stability of 
the bright and vortex solitons shown in figure [2] under small breathing oscillation started 
by a sudden change in scattering length a. This can be implemented experimentally by 
the FR technique. The initial soliton was created by imaginary-time propagation and 
subsequent dynamics generated by real-time propagation. In figure [3] we plot rms sizes 
(z) and (p) versus time. For the breathing oscillation of (p), we also plot a model result 
calculated from the theoretical frequency of twice the trap frequency uj± for an extreme 
cigar-shaped condensate using a hydrodynamic description [30] . For the bright soliton, 
the angular frequencies for (p) are uj p = 2uj± (num, var, and theory). For the vortex 
soliton, uj p = 1.94o;_l (num) and 2uj± (theory). The angular frequency for (z), controlled 
by the nonlinearity and not by any external trap, is very small in both cases. 

Now we investigate the collision between two bright and two vortex solitons of 
figure [2j Two such solitons are placed at z = ±25.6 at t — and are then advanced 
by real-time propagation of ([T]) with TV = 2000. Each soliton is attributed a velocity 
of about 5 mm/s towards center z = by including a phase factor exp(ivz) with v 
a constant in the initial wave function. The real-time simulation is terminated when 
the solitons reach z = ±25.6 at t = 11.22 ms. The collision dynamics is illustrated 
in figure [4] for bright and vortex solitons, where we show snapshots of contour plots of 
density 0, z ) t)\ 2 at different times before (0, 1.87 ms, 3.74 ms), during (5.61 ms) 
and after (7.48 ms, 9.35 ms, 11.22 ms) collision. The solitons come towards each other, 
interact at z = and then separate and come out practically unchanged. The quasi- 
elastic collision at such a large relative velocity of 1 cm/s demonstrates the solitonic 
nature and robustness of the bright and vortex solitons. Video clips of the collision 
dynamics of figure [4] are also prepared for bright and vortex solitons and contained 
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Figure 3. rms sizes (z) (lower panels) and (p) (upper panels) versus time during 
breathing oscillation of the bright (left panels) and vortex (right panels) solitons of 
figure [2] started by jumping a from 0.5 to 0.5025 nm and from 4 to 4.02 nm, respectively: 
numerical (num), variational [30] (var), theoretical (theory). 
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Figure 4. Contour plot of density \4>(x, 0, z, t)\ 2 of two colliding (a) bright and 
(b) vortex solitons of figure [2] with relative velocity 1 cm/s, before, during and after 
collision. 




Figure 5. Same as figures [2] (a) and (b) at time t = 11.22 ms after the collision 
considered in figure [4j 



in supplementary clips SI (also at http://www.youtube.com/watch?v=ds2CzZqN9M4) 



and S2 (also at http: / /www.youtube.com/watch?v=aQgLLhW08mE|), respectively. 



To illustrate the quasi-elastic nature of collision, we show in figures [5] (a) and (b) 
the 3D contour density plot of the final states at time t = 11.22 ms after the collision 
illustrated in figure [| for bright and vortex solitons, respectively. The similarity of plots 
in figures [5] (a) and (b) for the final states after collision with the initial states in figures 
2] (a) and (b) demonstrates the elastic nature of the collision. Specially, in figures [2] 
(b) and [5] (b) the radius of the vortex cores are practically the same determined by 
the fundamental healing length of the DBEC. To quantitatively compare the solitons 
before (t = 0) and after (t = 11.22 ms) collision, we show in Table I the numerical and 
variational results of chemical potential, rms sizes, and number of atoms of the initial 
and final solitons, in good agreement with each other. There is no exchange of atoms 
between solitons during collision. 

The quasi-elastic collision at high velocity is insensitive to the initial phase of the 
solitons. At low velocity, the interaction is sensitive to the initial phase. Two solitons 
placed side-by-side at rest along the z direction of opposite phase repeal, whereas of the 
same phase attract and form a stable bound soliton molecule [3]] . 
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Table 1. Numerical (n) and variational (v) /i, (z), (p), and N of the solitons of figures 
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Figure 6. Dynamics of soliton pair formation from contour plot of ID density versus 
time of (a) two bright and (b) vortex solitons of figure [2] placed side by side at rest. 




Figure 7. The initial 3D contour plots of the two (a) bright and (b) vortex solitons 
of Fig. [6] at t = and of the two soliton molecules formed of the two (c) bright and 
(d) vortex solitons at t = 400 ms. The vortex cores can clearly be seen in (b) and (d). 
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The dynamics of the formation of a soliton molecule from two equal-phase solitons 
of figure [2] is shown in figure |6| where we exhibit the contour plot of ID density 
d(z,t) = 2n f prfp|0(r, t)\ 2 versus time t and z. The initial positions of the bright 
(vortex) solitons are ±17.6 fim (±22.4 /xm). Due to dipolar attraction the solitons come 
close, coalesce and oscillate forming a stable bound soliton molecule [31j. In the case 
of the two bright solitons in Fig. [6] (a), after formation, the soliton molecule breaks 
into two solitons, then forms a molecule, thus oscillating between a molecule and two- 
soliton phases. In case of the collision of two vortex solitons, the final soliton molecule 
is a vortex molecule with zero density along the axial direction. This is illustrated in 
Fig. [7] where we show the 3D contour plots of the two interacting bright and vortex 
solitons of Fig. [6] at t — and 400 ms. A pair of bright solitons, bound together 
by a dark soliton, has been observed in dispersion-managed fiber optics [32]. Now it 
seems possible to observe the present soliton molecules in DBEC. It should be possible 
to consider an antivortex [33] in a DBEC using the present model based on a 3D GP 
equation. A problem of future interest is to consider the collision between a vortex and 
antivortex at low velocities. One may expect that a vortex- antivortex pair will annihilate 
the angular momentum quantum number, thus forming after collision a soliton without 
angular momentum, instead of forming vortex molecule of the type shown in Fig. [7] 
(d) with a hollow along the axial direction. Video clips of the molecule formation 
dynamics of figure [6] are also prepared for bright and vortex solitons and contained in 
supplementary clips S3 (also at |http: / /www.youtubexom/watch?v=tnnFOnbilgs| ) and 
S4 (also at |http:/ /www.youtube.com/watch?v=01AWIBzbrDs[ ), respectively. 

To conclude, we suggested the possibility of bright and vortex solitons in DBEC 
for repulsive short-range interaction for a < add and studied the statics (rms sizes and 
chemical potential) and dynamics (collision and breathing oscillation) of such solitons. 
Numerical results of stability phase plot and breathing oscillation of the bright solitons 
are in good agreement with variational calculations. The collision between two solitons 
is quasi-elastic at high velocities (~ 1 cm/s) and independent of initial phase. At very 
low (~ 0) initial velocity, two equal-phase solitons placed side-by-side at rest form a 
bound soliton molecule [3]] due to the dipolar attraction. With present technology, 
these solitons could be created in laboratory and the present theoretical results verified. 
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